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On the possibility of measuring the post-Newtonian 
gravitoelectric correction to the orbital period of a test 
body in a Solar System scenario 
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ABSTRACT 

The possibility of measuring the post-Newtonian gravitoelectric correction to the or- 
bital period of a test particle freely orbiting a spherically symmetric mass in the Solar 
System is analyzed. It should be possible, in principle, to detect it for Mercury at 
a precision level of 10 -4 . This level is mainly set by the unavoidable systematic er- 
rors due to the mismodelling in the Keplerian period which could not be reduced by 
accumulating a large number of orbital revolutions. Future missions like Messenger 
and BepiColombo should allow to improve it by increasing our knowledge of the Mer- 
cury's orbital parameters. The observational accuracy is estimated to be 10~ 4 from 
the knowledge of the International Celestial Reference Frame (ICRF) axes. It could 
be improved by observing as many planetary transits as possible. It is not possible to 
measure such an effect in the gravitational field of the Earth by analyzing the motion 
of artificial satellites or the Moon because of the unavoidable systematic errors related 
to the uncertainties in the Keplerian periods. In the case of some recently discovered 
exoplanets the problems come from the observational errors which are larger than the 
relativistic effect. 
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1 THE POST-NEWTONIAN 

GRAVITOELECTRIC CORRECTION TO 
THE ORBITAL PERIOD 

The geodesic motion of a test particle in the gravitational 
field of a spherically symmetric body of mass M is con- 
sidered. In the framework of the Einstein's General The- 
ory of Relativity the post-Newtonian gravitoelectric accel- 
eration of order Q(c~ 2 ) experienced by the test particle is 
Jjoos fc Grafarendlll99ll) 
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where G is the Newtonian gravitational constant, c is the 
speed of light in vacuum, r and v are the position and ve- 
locity vectors, respectively, of the test particle, r is the stan- 
dard isotropic radial coordinate , not to be confused with 
the Schwarzschild radial coordinate r = r[l + GM/ (2c 2 r)] 2 . 



* E-mail: lorenzo.iorio@libero.it (AVR) 

1 It is the coordinate used in t he force models of the post- 
Newtonian equations of motions <EstabroolJll97ll) adopted for 
the computations of the planetary ephemerides by, e.g., the Jet 
Propulsion Laboratory (JPL). 



The gravitoelectric acceleration of eq.Q is the same which 
is res ponsible for th e Einstein secular advance of the peri- 
helia (lEinstemlll915l) of the Solar S ystem planets which are 
known at a 10 -4 level of accuracy jPitievall2001allcl . 

We are interested in the consequences of eq.Q on the 
orbital period of the test particle in order to see if they can 
be detected in suitably designed experiments in the Solar 
System arena. Let us first consider circular orbits for which 
r ■ v — and r = a, where a is the semimajor axis of 
the orbit; indeed, the orbits of the planets and of many of 
the best tracked Earth artificial satellites have rather small 
eccentricities e. Subsequent sensitivity analyses will show 
the extent to which this approximation can be considered 
adequate. In this case the acceleration given by eq.Q is 
directed radially and the radial equation of motion reduces 
to 
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With v — ad(j>/dt, the time required for the azimuthal angle 
<J) to pass from to 2ir becomes 
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where P^ - 1 is the unperturbed Keplerian period Note that 
cq.JSJ is consistent with the result of ( Mashhoon et all 
2001). Contrary to the g ravitomagnetic correct ion to the 
Keplerian orbital period l|MasIihoon et al] l20011. which is 
generated by the proper angular momentum J of the cen- 
tral mass, p' GE ) is insensitive to the direction of motion of 
the test particle around its orbit and to its inclination with 
respect to the equatorial plane of the central body. 



2 A TENTATIVE ERROR BUDGET 

Is it possible to measure p( GE ' by analyzing the motion of 
the planets around the Sun? Although their eccentricities 
are rather small being of the order of 10 -2 — 10 -3 , apart 
from Mercury for which e = 0.2056, it turns out that, as we 
will show below, the available observational accuracy does 
not allow to neglect the impact of the non-circularity of their 
orbits. For e ^0 the exact e xpression of the post-Ne wtonian 
gravitoelectric correction is ijMashhoon et alJ^OOlf ) 
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where /o is the true anomaly at the initial epoch of the 
planet. The factor 
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runs from 1.65 for f = 0° to -0.10 for f = 180°; in the 
case of Venus, which has e = 0.00677, it changes from 1.02 
to 0.97. In the following, in order to sketch a sensitivity 
analysis, we will work with the initial epoch for which the 
factor of eq.Q is equal to 1, i.e. we will assume eq.(JSJ. 



2.1 The observational sensitivity 

The azimuthal angle <j> is the usual right ascension a of 
a spherical coordinate system in an astronomical reference 
frame. 

We will refer to the International Ce lestial Reference 
Frame 2 (ICRF) jMcCarthv fc Petit) 12003) . Moreover, the 
right ascension is one of the direct observables in planetary 
motions. Since we are interested in the times when the right 
ascension of a planet crosses the {x,z} plane of ICRF, i.e. 
when it is zero, it is of the utmost importance to see if the 
present uncertainty in the stability of the ICRF's axes would 
allow for a detection of the gravitoelectric correction or if it 
is too large and would mask it. We can reasonably assume 
that the right ascension of a planet advances under the effect 
of the post-Newtonian gravitoelectric force over an angular 
interval of a = 2n, by a quantity 5a denned by 

P< GE > 3GM 

(6) 
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where a is the semimajor axis. From eq.© it is appar- 
ent that the inner planets have to be considered in or- 
der to maximize the effect. For Mercury, whose semimajor 

2 It is the re alization of the Int ernational Celestial Reference Sys- 
tem (ICRS) lArias e t al.ll99 5|) bv means of the estimate s of a set 
of extragalactic sources <Ma et aDll.99St iFev et aljEo04l) . 



axis amounts to a = 0.38709893 A.U., P (GE) = 0.29 s (if 
the numerical factor of eq.JKJ is assumed equal to 1) and 
P (0) = 7.60055184 x 10 6 s; eq.@ yields 5a/ a = 3.8 x 10 -8 . 
Then, 5a = 50 milliarcseconds (ma s in the following). Ac- 
cording to llMcCarthv fc PetTrll2004h . the uncertainty in the 
ICRF axes amounts to a a = 0.02 mas. This yields an accu- 
racy of 4 x 10 -4 . Then it appears clear that with such an 
accuracy 3 the circular case approximation cannot be consid- 
ered adequate. 

Such evaluations hold for a single orbital revolution 
only; a measurement like that proposed here must be per- 
formed over a large number of planetary transits. This would 
greatly increase the observational accuracy. 

Our knowledge of the orbital motion of Mercury will 
improve thanks to the future missions Messenger 4 , which 
has been launched in the summer 2004 and whose encounter 
with Mercury is scheduled for 2011, and, especially 5 , Bepi- 
Colombo 6 , which is scheduled to fly in 2010-2012. A com- 
plete error analysis for the range and rang e-rate measure- 
ment s to BepiColombo can be found in iless fc Boscaeli 
l200lT) . According to them, a two orders of magnitude im- 
provement in the Earth-Mercury range, which is accurate 
to hundreds of meters now, should be possible. Accord- 
ing to a more conservative evaluation by E.M Standish 
(JPL)(Standish, private communication 2004), improve- 
ments in the Mercury's orbital parameters might amount 
to one order of magnitude. In regard to the proposed test, 
the Mercury transits which could benefit from Messenger 
and BepiColombo will be of the order of 4, since both the 
spacecraft should nominally orbit Mercury for 1 year and the 
orbital period of Mercury around the Sun is almost equal to 
88 days. 



2.2 Some systematic errors 

2.2.1 The errors in the Keplerian period 

The value of the Keplerian period, evaluated from the es- 
timated semimajor axis a and Sun's GM, must be sub- 
tracted from the data record in order to single out the 
post-Newtonian effect (and the other classical and post- 
Newtonian perturbations, of course, which, in this case, 
would represent the noise). Then, let us see if the system- 
atic errors in the Keplerian period are smaller than the post- 
Newtonian gravitoelectric correction which we are interested 
in. This is a very important point because, opposite to the 



3 It maybe interesting also to note that such an accuracy would 
not allow for a detection of the gravitomagnetic correction to 
the Keplerian period 2ttJq/c 2 M which, in the case of the Sun, 
amounts to 4x 10 -6 s by assuming (Piipcrs 2003) Jq = 1.9 X 10 41 
kg m 2 s _1 . 

4 See on the WEB |http://messenger.jhuapl.edu/| and 
|http: //discovery.nasa.gov/msssenger.html| 

D While the spacecraft trajectory will be determined from the 
range -rate data, the plan et's orbit will be retrieved from the range 
data iMilani et alJl2002l) . In particular, the determination of the 
planetary centre of mass is important to this goal which can be 
better reached by a not too elliptical spacecraft's orbit. The rel- 
atively moderate ellipticity of the planned 400 X 1500 km polar 
orbit of BepiColombo main orbiter, opposite to the much more 
elliptical path of Messenger, is, then, well adequate. 

6 See on the WEB htt p://sci.esa. int/science-e/www/area/index.cfm?fareaid=30| 
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observational error, the impact of this source of systematic 
bias could not be reduced by observing many orbital revo- 
lutions. We have 
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By assuming 7 <t m^c = 1.87 m .PitievalE-OOlaT) and gg m = 
8 x 10 9 m 3 s -2 iStandishlll99Eh . we have cr pm ^ 3 x 10" 4 
s + 2 x 10~ 4 s. The NASA Messenger and the future ESA 
BepiColombo missions should allow to further reduce cr p (o) 
because it will yield a better knowledge of the parameters 
of the Sun-Mercury system. 

The effect of the quadrupole mass moment J2Q of the 
Sun would induce a correction 
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By assuming the rang e J20 = (2 ± 0.4) x 10 7 
.Pireaux fc Rozelotlliooai . it amounts to (6 ± 1.2) x 10~ 4 
s. Note this error should be reduced to ~ 10~ J s if J20 will 
be measured with an accuracy of the order of 10 -9 by Bepi- 
Colombo. These sources of constant bias, which are at the 
same level of the one-revolution observational error previ- 
ously outlined, set the limit of the obtainable accuracy over 
a record spanning over many planetary revolutions. 



2.2.2 The direct planetary secular perturbations 

Another source of perturbations on the Mercury's right as- 
cension is represented by the gravitational perturbations in- 
duced by the other major bodies of the Solar System. Let 
us calculate the secular effects induced by some of the other 
planets. The perturb ative effect of the planet m 7 - on the 
planet riii is given by llBoccaletti fc PucaccollT999l) 
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It turns out that the second term in eq.@ does not in- 
duces secular perturbations. After expressing the first term 
of eq.© in terms of the orbital elements of the i-th and 
j-th bodies and averaging it over one period of the mean 
longitudes A* and X J it can be obtained, for the largest con- 
tribution which does not contain the terms of second order 
in the eccentricities and the inclinations 
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The nominal values of eo. ljKty are larger than the gravito- 
electric effect: the major contributions come from Jupiter 
(-5.97 s), Venus (-5.70 s), the Earth (-2.64 s) and Saturn 
(-0.29 s). If the estimates of eg. 1101 are correct, this would 
not pose problems because the systematic errors induced by 
o~ a i and cr Gm j would be far smaller than the gravitoelec- 
tric effect. Inde ed, by assuming for Jupiter cr a j up = 675.1 m 
JPitievall200Tah and 8 <7 Gm .iu P = 2 x 10 9 m 3 s -2 .Jacobsonl 



7 In fact, the formal, statistical error is 0.187 m: in this case 



Jp(0) 



< 3 X 10 5 s. However, the real error bound in a might be 



one order of magnitude larger. Note also that according to some 

authors the solar GM should be considered as fixed. 

8 See on the WEB http://ssd.jpl.nasa.gov/sat_gravity.html 



2003) the uncertainty in its contribution would amount to 
10 -7 s only. The nominal values of the corrections induced 
by Jupiter allow to guess that the impact of ma ny of the 
aster oids, whose masses are not accurately known .Standisbl 
2002), should be negligible: indeed, from eg. (HOP we can as- 
sume that p( astcroid ) /p(Jup) 



, (asteroid) / m (J u P) 



2.2.3 The indirect errors in the Keplerian period due to 
the short-term variations of the semimajor axis 

A possible source of aliasing for Mercury could be repre- 
sented by the indirect perturbations 
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induced by short-periodic effects on a, i.e. those effects 
which are not averaged over one orbital revolution of the 
various planets. They are 

• The indirect effects on P' ' induced by the high- 
frequency perturbations on Mercury's semimajor axis due to 
post-Newtonian gravity itself. They can be calculated from 
eq.jy and the Gauss perturbative equation for the rate of 
the semimajor axis. It turns out that they are 

A « (GE) = ^f)2N ( COs/0--COs/) + 

+ 10e(cos 2 / - cos 2 /o)] +0(e 3 ). (12) 

For Mercury their nominal amplitudes are of the order of 
4 km; by assuming <tgm = 8 x 10 9 m 3 s -2 the error in 
them is of the order of 10 -7 m. We can, then, conclude that 
the indirect effects due to the post-Newtonian short-period 
shifts in the semimajor axis are negligible. 

• The indirect effects on P' ' induced by the high- 
frequency perturbations on Mercury's semimajor axis due 
to the classical N-body perturbations. In order to get just 
some orders of magnitude, let us see what could be the per- 
turbation induced on the semimajor axis of Mercury by a 
planet whose orbital elements will be marked with . It turns 
out that, for the term of lowest degree \l\ + \l | + \m\ + \m 
which must be ^ 2 because in the Lagrange planetary equa- 
tions only the first partial derivatives appear, we can write 
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x sin(jA + j A + Izu + I w + mfl + m Q )dt, (13) 

where j^O, zu = cu + is the longitude of perihelion 
and is the longitude of the ascending node. The condition 
j + j +1 + 1 + m + m = must be, in general, fulfilled 
together with that which states that m + m must be an even 
number (positive, negative or zero). Since now we are not 
looking at the secular terms, the supplementary condition 
l + l +m + m 7^ must be fulfilled as well. It turns out that 
mixed terms proportional to, say, ee 2 contribute to eg. 11311 . 
For Jupiter, the factor 2GM /n 2 aa amounts nominally to 
8.2258248 x 10 6 m. The uncertainty due to a Gm ju P = 2 x 
10 9 m 3 s" 2 is of the order of 10" 1 m, while a a ju P = 675.1 
m yields an error of the order of 10 -3 m. It seems, then, 
reasonable to conclude that the indirect perturbations on 
the semimajor axis should not pose severe limitations to the 
proposed measurement. 
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3 THE (IM)POSSIBILITY OF A 
MEASUREMENT WITH EARTH 
ARTIFICIAL SATELLITES AND THE MOON 

Let us now see if the gravitoelectric correction p' GE > could 
be measured by analyzing the motion of Earth artificial and 
natural satellites whose eccentricities are rather small. 

From an observational point of view, by repeating the 
reasonings of Section ^. ll with eq.© and by considering that 
the accuracy in the orientation of the axes of the Interna- 
tional Terrestrial Reference System 9 (ITRS) is of the order 
10 of 3 mas, it would be possible to measure p( GE ' over a 
sufficiently high number of orbital revolutions. Indeed, for a 
typical satellite orbit with a — 7 x 10 6 m the post-Newtonian 
shift over 2-zr would amount to 1 mas. 

However, by assuming, e.g., a = 1.2270 x 10 7 m, as for 
LAGEOS whose rms post-fit residuals amount to a a ~ 10 -2 
m or better, the systematic error on the Keplerian period 
P (0) amounts to lxlO" 5 s, while P (GB) is of the order of 
7 x 1CP 6 s only. It is important to note that this uncertainty, 
opposite to the observational one, could not be reduced by 
observing a large number of transits: indeed, the 1-cm level 
of accuracy in estimating the semimajor axis comes just from 
a data reduction over a time span of many orbits. 

Let us, now, consider the possibility of using the very ac- 
curate present and future data from the Lunar Laser Rang 
ing (LLR) technique for the Moon 



Recalling that a 



3.84400 x 10 s m and e R 



= 0.0554 



(Standish 2001), the approximation of circular orbit would 
be adequate and p( GE ' = 4 x 10 -5 s almost independently 
of /o; the gravitoelectric advance over 2n amounts to 7.5 
mas. By assuming the present day centimeter accuracy in 
knowing the lunar orbit, the systematic error related to the 
Keplerian period, which cannot be reduced by accumulating 
a large number of revolutions, would amount to 1 x 10~ 4 s. If 
the millimeter accuracy level will be reached l|Williams et alJ 
l2004h . the systematic error will be reduced down to 1 x 10 
s, i.e. 25% of the gravitoelectric correction. 

The major limiting factor both for the Moon and the 
artificial satellites is the impact of the uncertainty in the 
Earth's GM which amounts to <jgm^ — 8 x 10 5 m 3 s -2 
(McCarth y fc Petitil2004f l. Indeed, the ratio of the error in 
the Keplerian period due to <jgm 9 to the gravitoelectric 
term amounts to 



r p(°)lGM ffi C 2 (aGM @ ) 
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a = (1.5084 x 10" 
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Then, for a typical satellite orbit with a = 7 x 10 6 m this 
systematic bias would be close to 100%. For the Moon it 
would be 5800%. 



9 Its re alization, the International Terrestrial Reference Frame 
(ITRF) iMcCarthv fc PetiJ20()l) . is obtained by estimates of the 
coordinates and velocities of a set of observing stations on the 
Earth. 



4 THE (IM)POSSIBILITY OF A 

MEASUREMENT WITH EXOPLANETS 

A potentially interesting category of celestial objects which 
fit the conditions of validity of eq.© is represented by 
many exoplanets 11 recently discovered with the transit 
method. Indeed, they are giant planets following circular 
trajectories with orbital periods of the order of 1 day, i.e. 
very close to their star. On the other hand, the accu- 
racy of the measurement of their periods with photomet- 
ric techniques is continuously increasing. L et us look at 
the very recently analyzed OGLE- TR- 132b feouchv et alJ 
120041: iMoutou et aTll2004D . The relevant orbital parame- 
ters are M/Mq — 1.35 ± 0.06 (adopted for the star), 
a/A.U. = 0.0306 ± 0.0008, m (c) /m (Jup) = 1.19 ± 0.13 
and P (obs) = 1.46003645 x 10 5 s ± 0.518 s (measured for 

the planet over many orbital revolutions). Unfortunately, 
p(GE) = Q95 g Qnly 

The situation is even worse in the case of the pulsars' 
planetary systems. For, e.g., the 0.020 m® planet orbiting 12 
PSR 1257+12 at 0.19 A.U. the uncertainty in the observed 
orbital period is 259.2 s while P (GE) is 0.24 s. 



5 CONCLUSIONS 

In this paper we have conducted a preliminary sensitiv- 
ity analysis about the possibility of measuring the post- 
Newtonian gravitoelectric correction to the orbital period 
of a test body in a Solar System scenario. The conclusions 
are summarized in Table 1. 

The post-Newtonian shift of Mercury's right ascension 
would amount to 50 mas over one orbital revolution. Its 
orbit has been assumed circular just in order to perform 
a preliminary sensitivity analysis. In fact, the obtainable 
observational accuracy does not allow to neglect the impact 
of the factor of eq.JSJ which depends on the eccentricity and 
the true anomaly at the initial epoch. 

Such an effect should be measurable, in principle, for 
Mercury with a 10~ 4 accuracy. This level is mainly set by 
the systematic errors induced by the mismodelling in the 
Keplerian period which should be subtracted. Instead, the 
precision with which the axes of the International Celestial 
Reference Frame are known (0.02 mas) can be assumed as 
representative of the observational accuracy. It yields a 10~ 4 
accuracy over one orbital revolution only. Of course, it could 
be increased by observing as many transits as possible. Fu- 
ture improvements in astrometry and in the knowledge of 
Mercury's orbital parameters from the Messenger and Bepi- 
Colombo mission should allow to further increase the pre- 
cision of this test by reducing the systematic errors. Bepi- 
Colombo is also expected to measure the Solar quadrupole 
mass moment Ji with an accuracy of few parts in 10~ 9 . 

The present-day level of accuracy in knowing the or- 
bits of the best accurately tracked Earth artificial satellites 
rules out the possibility of measuring such post-Newtonian 
effect in the terrestrial gravitational field. For the Moon, a 
millimeter accuracy in knowing its orbit, a goal which could 
be reached in the next future, would yield a 25% systematic 

11 See on the WEB http://cfa- www.harvard.edu/planets/catalog.html 



See on the Internet http://hpiers.obspm.fr/iers/bul/bulb/explanator;f.4it$Efc http://www.obspm.fr/encycl/1257+12.html 
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Table 1. Observational and systematic errors affecting the measurement of the post-Newtonian gravitoclcc- 
tric correction p( GE ) to the orbital period of Mercury, the Moon, the Earth artificial satellite LAGEOS, 
an exoplanet of the system OGLE-TR-132b and a planet of the PSR 1217+12 pulsar. All figures are in s. 
The circular orbit approximation has been adopted. The observational accuracy is to be intended over one 
orbital revolution. It is retrieved from the ratio of the accuracy, in mas, with which the axes of the ICRF 
and ITRS are known to the gravitoelectric shift over one orbital revolution. Contrary to the systematic er- 
rors, the observational errors can be reduced by observing a sufficiently large number of orbital revolutions. 
Note that the observational uncertainty for the exoplanet refers to many orbital revolutions. The systematic 
errors are due to the Keplerian period p(°) , the indirect effect AP' 8 ' on it due to the perturbations in the 
scmimajor axis, the perturbation induced by the quadrupolar mass moment of the central source P( j a), the 
direct perturbation induced by the secular N-body interaction p(p lanets ) . T n regard to the indirect effects on 
the Keplerian period, f AP (o) |ge and <T Ap ( ) Ipianets are induced by the short-period post-Newtonian and 
Newtonian effects on the scmimajor axis. 





Mercury 


Moon 


LAGEOS 


OGLE-TR-132b 


PSR 1217+12 


p(GE) 


2.9 x 1CT 1 


4 x KT 5 


7 x itr 6 


9.5 x KT 2 


2.4 x itr 1 


fobs 


1 x 10" 4 


1.6 x icr 5 


2 x itr 5 


5.18 x io -1 


2.592 x 10 2 


O-p(O) |ct 

">(<)) Igm 

<7p(planets) 

°"ap<°) Ige 

O'AP(O) Ipianets 


3 x itr 5 
2 x itr 4 
1.2 x icr 4 

< itr 7 
<io- n 

< icr 7 


1 x icr 4 

2 x icr 3 


i x itr 5 
l x icr 5 







bias due to this source of error. Instead, the present-day cen- 
timeter level in LLR would not allow to sufficiently reduce 
the systematic error due to the uncertainty in the Keplerian 
period. However, the uncertainty in the present-day knowl- 
edge of the geocentric gravitational constant GM ffl would 
induce a systematic error of more than 100% for every ter- 
restrial artificial or natural satellites. 

The exoplanet scenario is not known with a sufficiently 
accuracy from the point of view of the observational errors 
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